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' Abstract 

■ By combining results of Mizumachi on the stability of solitons for the Toda lattice with a 

simple rescaling and a careful control of the KdV limit we give a simple proof that small ampli- 
tude, long-wavelength solitary waves in the Fermi-Pasta-Ulam (FPU) model are linearly stable 
and hence by the results of Friesecke and Pego that they are also nonlinearly, asymptotically 
stable. 
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1 Introduction 

In a series of four recent papers Friesecke and Pego ([T], [2], [3], [4]) made a detailed study of the 
existence and stability of solitary wave solutions of the Fermi-Pasta-Ulam (FPU) system: 



q,^V'{q, + ^-q,)-r{q,-q,^^) jeZ (1) 



which models an infinite chain of anharmonic oscillators with nearest-neighbor interaction potential 
V. If we make the change of variables Vj = qj+i ~ qj and pj = qj, the state variable u = (r,p) 
satisfies a system of first order Hamiltonian ODEs, 

oo ■ 

O ■ Uf = JH'{u) (2) 



where the Hamiltonian H is given by 

H{r,p)=Y.\pl + ^^''k) (3) 

The symplectic operator J is given by J = ^ ^ ^g-i q ^ ^ where S* is the left shift on bi-infinite 

sequences, i.e. {Sx)n '■= Xn+i- The problem is well posed in each P space, but for concreteness and 
simplicity we work in l'^. Throughout the paper we shall assume that the interaction potential V 
satisfies the following 

yeC*; y(0) = y'(0) = 0; X^"(0)>0; ^"'(0) 7^ (4) 

In the first paper in the series, [1], Friesecke and Pego prove that the system ([2]) has a family of 
solitary wave solutions which in the small amplitude, long-wavelength limit have a profile close to 
that of the KdV soliton. More precisely they show: 

Theorem 1.1 (Theorem 1.1 (b) from [T], restated). Assume that V satisfies (jj]) and that c > 
\/v'''JO) is sufficiently close to ^V"(0). Then there exists a solution to the wave profile equation 
for FPU: 

cr'^{x) ^{S + S-^ - 2/)V'(re) re(±oo) = 



1 



which in addition satisfies 

1 



7?-^ - 



< Ce' (5) 
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where 4)i{x) := ^„/^°^^ (^^sech{^x))^ is the KdV soliton and e := y^^^ 

Remark 1.2. In fact, in it is only proven that the traveling wave profile is close to the (rescaled) 
KdV soliton in the norm. The strengthening of the estimate to hold in the norm was done 
in ig/. 

Remark 1.3. The first general results about the existence of traveling waves in these general FPU 
type systems were obtained by Friesecke and Wattis via variational methods. 

In the second paper in this series, [2], the authors use the method of modulation equations to 
prove that if the sohtary waves are hnearly stable, they are nonlinearly stable as well. More precisely, 
suppose that the following linear stability condition is satisfied: 
Define 



-1 



^((?-,p),(/3,7r)) ( Pk+jPj+ X! ^k+jT^j] (6) 

and define 



\k— — oo k- 



-oo 



||x||^=^e^'^^x|. (7) 
fcez 

Hypothesis L: There are positive constants K and [3' and cq > \/V"(0) such that whenever \/V"(0) < 
c* < Co, and w is a solution of the linear equation 

dtw = JH"{rLc,{- - cj))w (8) 

with ||u'(to)||a < oo and such that 

i^idzUc,{z)\z=--c,t,w{to)) = Uj{dcUc{z)\c=c,,z= — c,t,w{to)) = (9) 

holds, then the estimate 

e--'*\\w{t)\U < i^e-'^'(*-^)e— *^||^(s)|U (10) 

holds for all t > s > to. 

Theorem 1.4 (Theorem I.f from [5], restated). Assume that V satisfies Assume that the 

semigroup generated by the FPU model linearized about its solitary wave satisfies Hypothesis L. Then 
the solitary wave Uc is stable in the nonlinear system ([5]) in the following sense: Let (3 G (0,/3'). 
Then there are positive constants Co and 5q such that if for some 6 < 5q and 7* S M the initial data 
satisfy 

ho-"c.(--7*)ll < and \\e<-'<'\uo - u,,{- - < 5 

then there is an unique asmyptotic wave speed Coo o,nd phase 700 such that 

|coo - c*| + |7oo - 7*1 < Co5 

and 

\\u{t,-)~Uc^{- - Coot-7oo)|| < CoV^ t > 0, 

and in addition 

||ga(.-c^t-7oo) .) _ ^^^(. _ _ ^^)) II < c^se-f* 
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The last two papers in the series, [3], [1], are devoted to verifying that the hnear estimate ([TU|) 
holds for the solitary waves constructed in [T] . 

In [3], Friesecke and Pego construct a type of Floquet theory to prove estimates like (fTO)) . The 
reason that one needs such an approach is that because of the discreteness introduced by the lattice, 
the linearized equation ([5]) is not autonomous in a frame of reference moving with the traveling 
waves, but only periodic (with a spatial translation.) Finally, in paper [1], Friesecke and Pego verify 
that the solitary waves constructed in [T] satisfy the hypotheses of their Floquet theory and hence, by 
the results of [2] are asymptotically stable. This last step involves, among other things, the fact that 
these solitary waves are well approximated by the KdV soliton and the fact that the linearization of 
the KdV equation about its solitary wave is well understood. 

In this note, we give a simple alternative proof of the estimate pO)) which avoids the use of the 
Floquet theory and spectral analysis of [3] and [4|. 

Our proof combines three observations: 

1. The linear stability of the soliton for the Toda lattice, established by Mizumachi and Pego by 
the construction of an explicit Backlund transformation in [7]. 

2. A transformation of the original FPU equation ([2]) into a form in which we can prove that its 
solitary wave solution is close to that of the Toda lattice. 

3. A careful control of the way in which various quantities depend on the small parameter e. 

We note that the last two points were originally developed in our study of counter-propagating 
2-soliton solutions of the FPU model [5] . 

We now explain these three points in more detail. Recall first that the Toda lattice is the special 
case of the FPU model with potential function 

V{x) ^ {e-^^ + X-1) . (11) 

(Throughout this paper, quantities with tildes will refer to the Toda model.) Note that V satisfies the 
hypotheses of [1] so the results of that paper imply that the Toda model has a family of solitary waves 
close to those of the KdV equation. Of course in the case of the Toda model these solutions were 
explicitly constructed by Toda in the 1960's and indeed the Toda model is one of the classic examples 
of a completely integrable, infinite-dimensional Hamiltonian system. The stability of the Toda soliton 
can also be analyzed in a very direct fashion. By constructing a Backlund transformation which 
conjugates the linearization of the Toda model about its soliton to the linearization of the Toda lattice 
about zero, Mizumachi and Pego proved that the linearized Toda equation satisfies Hypothesis L 
and hence, by the results of [5], that the Toda soliton is asymptotically stable. In [S] we extended 
this result by showing that the constant K in PT7)1 can be chosen uniformly in c. 

Remark 1.5. Although we will be most interested in these results in the long-wavelength, small 
amplitude regime studied by Friesecke and Pego, in fact the results of apply to Toda lattice 
solitons of arbitrary size. 

In order to extend the estimate on the linear decay from the semigroup of the linearized Toda 
equation to the linearization of the general FPU model satisfying hypothesis (HI), we first make 
use of the following simple: 

Lemma 1.6. Without loss of generality we may assume that the potential energy function V in 
( or equivalently ([3]) ) satisfies 

v"io) = 1 , y"'(o) = -1 . 



Proof. To sec this simply note that if q{j,t) is a sohition of ([T|) with a potential energy function 
V that satisfies hypothesis (HI), then V can be written as V{x) ~ ^ax^ + + 0{x'^) for some 
a > and 6^0. If we now define q{j-t) = aq{j, (3t), then q solves the FPU equations with potential 
function V{x) = ^ax^ + |6(/3^/a)a;'^ + 0{x'^), so choosing a and (3 appropriately insures that the 
lemma holds. □ 
From now on we will assume that the potential function V satisfies this normalization. Note 
that with this normalization the potential 1/ in ([1]) differs from the Toda potential only by terms of 
0{x'^) or higher. With this observation the existence results of [1] imply: 

Proposition 1.7. Let Uc by the profile of the solitary wave of the FPU model ((!]) with speed c, and 
let Uc he the profile for the special case of the Toda potential. Define = dzUc{z) and ^2 = dcudz) 
and let ^12 be the corresponding quantities for the Toda lattice. Then there exists eo,C > such 
that for < e < eo, one has the estimates 

\\uc - UcWei < _ < (j^i (^12) 

Ui-ii\\il<Ce''^' Ui-^i\\i-<Ce' (13) 
116 - i2\\el < Ce^l^ 116 - UV^ < Ce^ (14) 

Remark 1.8. The half-powers of e that occur in the estimates of the £^ norms are a consequence 
of the scaling of the functions u^. 

Proof. The inequality ||mc — Mc||^=^ < Ct^ follows immediately from ([5]) because to leading (e^) order 
both Uc and Uc agree with the KdV soliton. The estimate on \\rLc — UcWt^ then follows from this 
estimate because of the prior remark about the scaling of the £^ norms. The estimates on 6 ~ ?i 
then follow from these two since Theorem 11.11 shows that a derivative of the solitary wave profile 
with respect to the spatial variable gains exactly one power of e. The estimates for 6 ^ 6 follow in 
a similar fashion. For more details see |S] . □ 
With these estimates in hand we consider the semi-group generated by 

dtv = JH"{uc)v = JH"{uc)v + J [h"{uc) - H"iuc)) v . (15) 



The idea is now to treat the term J \^H"{uc) — H"{uc)j w as a perturbation of the Toda semi-group. 
Recalling that H and H differ only at quartic order and that Uc and Uc are both of order ©(e^) and 
differ only by terms oi 0{e^) we have 

\\J (h"{uc)-H"{uc))v\\,2 < C {\\uc - UcWe^ + iWucWi^ + WucWt^f) \\v\\ei (16) 

The other fact we must deal with is that v E E'^ = {v \ v) = 1-^^(61 v) ~ 0}, while our decay 
estimates on the Toda semigroup hold only if the semigroup acts on vectors v E E'"^ ^ {v \ a;(6, v) = 
^{S,2tv) = 0}. To cope with this difference we define the projection operator 

^ / ^^(6,6)w(6,i') +^(6, 6)^(6,-^)1 r uj(ii,v) ~ 

Qv = v-\ 6 - ~ r 6 , (17) 

V ^(6,6)^ / ^(6,6) 

which maps i"^ to Eg. Using the estimates in Proposition 1 1.71 . we have: 

Proposition 1.9. [Lemma 4-4 JBl; simplified] There exists eo, C > such that ifO<e<eo, then 

\\Qv\\e2< C\\v\\e2 . 
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Furthermore, if v €z -B* , then 

Wil - Q)v\\r2 < Ce'/^v\\,. . 

Proof. The key factor in the proof of this proposition is that the estimates of Theorem 11.11 make it 
possible to evaluate the leading order in e behavior of the symplectic products uj{^j,^k) and uj{^j,^k)- 
This was first used in [T] and was utilized repeatedly in [5]. For instance, one has = 0, 

while tj(^2, $2) = C22e~^ + C(e~^) and the cross termw(^i,f2) = ci2e + C'(e^) with the constants C12 
and C22 both non-zero. Similarly, the leading order behavior in e of the norms of ^1 and ^2 can be 
computed by relating them to the derivitives of the profile of the KdV soliton using the estimates of 
Theorem ll.il The same estimates also hold for the symplectic inner products of ^1.2 and with these 
estimates the first bound in the proposition follows immediately. 

The second estimate follows by rewriting the projection operator as 

uj{i2,t.i)[I - Q)v = - , 4i-t^Ui,w)42 (18) 

V ^(6,6) / 



t^(6,Cl)^(l2 " 6,w) +^{l2,i2)^i.il,v) 



where the last step used the fact that since v G uj{^i,v) = uj{^2,v) = 0. But now note that 
each term on the right hand side contains a factor of either ^1 — ^1 or ^2 — £,2 and these are small 
due to the estimates in Proposition 1 1.71 . With the aid of these estimates the second estimate in the 
proposition follows in a straightforward fashion. For more details see the proof of Lemma 4.4 in [5] . 

□ 

Corollary 1.10. There exists eg > such that if < e < eg and v € i?*, then 

\\v\\,2<2\\Qv\\,2 . 
Now write v{t), the solution of p5|) as 

v{t) = S{t, O)w(O) + ^ Sit, s)J (h"{uc) - H"{uc)) v{s)ds , (19) 

where S is the evolution semigroup generated by the linearized Toda system. Then 

Qv{t) = S{t,0)Qv{0) + J S{t,s)Q(^j(^H"{u,)^H"{u,)'^yis)ds , (20) 

Taking the norm of both sides and using the estimates on the Toda semigroup (jlOp we find: 

||gt;|l,. < Ke-'^''\\vm,.+K J\-'^^''^'-^^\\q(^j(^H'\u,)-H'\u,^^^ (21) 

pt 

< Ke-^''''\\v{0)\\e2 +K2e^ e-'"'''-*-'^\v{s)\\pds (22) 

Jo 

< Ke-^''*\\v{0)\\e2 +K3e^ e-'"'^'-'^\\Qv{s)\\e2ds (23) 



^0 

Note that the constants K, K2, and are all independent of e. Now setting (pit) ~ supQ<^<( e} '^^'^||Q7j(t)|| 
for some < b' < b and taking the supremum in the above equation we have 

(/)(t) < K\\v{0)\\e2 + Kse^e-'^^-'''^''' f e-'^''-'''^''ms)ds (24) 
< K\\viO)\\e2+K3eq^{t) . (25) 
Thus, if e is sufficiently small, we conclude that (j){t) is uniformly bounded for all t and hence: 



Proposition 1.11. There exists K' > 0, independent of e, and eq > such that i/ < e < eq ^'^'^ 
v{t) €z i?" is a solution of psp then 

\\v\\e2<K'e-'''^''\\vm,. . 

This verifies tliat the hnearized FPU semigroup satisfies Hypothesis L and hence by the results 
of [2] that the FPU sohtary wave is asymptoticaUy stable and the small amplitude, long-wavelength 
regime. 

Remark 1.12. We note that in the proof of Provosition U. 1 1\ it is important to carefully control the 
dependence of the semi- group on the parameter e. It is not surprising that a perturbation arguments 
permits one to extend the results of Mizumachi and Pego to solitary waves in small perturbations 
of the Toda model. Indeed, this was already noted in . What we do find noteworthy is that this 
simple argument can cover all FPU solitary waves in the KdV regime, and this requires the detailed 
study of the small e asymptotics contained in Propositions | _? . 7| and \1.9\ . 

Acknowledgements: This work was funded in part by the National Science Foundation under 
grants DMS-0603589 and DMS-0405724. Any findings, conclusions, opinions, or recommendations 
are those of the authors, and do not necessarily reflect the views of the NSF. 

References 

[1] G. Friesecke and R. L. Pego. Solitary waves on FPU lattices. I. Qualitative properties, renor- 
malization and continuum limit. Nonlinearity, 12(6):1601-1627, 1999. 

[2] G. Friesecke and R. L. Pego. Solitary waves on FPU lattices. II. Linear implies nonlinear stability. 
Nonlineanty, 15(4):1343-1359, 2002. 

[3] G. Friesecke and R. L. Pego. Solitary waves on Fermi-Pasta-Ulam lattices. III. Howland-type 
Floquet theory. Nonlinearity, 17(l):207-227, 2004. 

[4] G. Friesecke and R. L. Pego. Solitary waves on Fermi-Pasta-Ulam lattices. IV. Proof of stability 
at low energy. Nonlinearity, 17(1):229-251, 2004. 

[5] Aaron Hoffman and C. Eugene Wayne. Counterpropagating two-soliton solutions in the FPU 
lattice. Preprint; arXiv:0806.1637, 2008. 

[6] Gero Friesecke and Jonathan A. D. Wattis. Existence theorem for solitary waves on lattices. 
Comm. Math. Phys., 161(2):391-418, 1994. 

[7] Tetsu Mizumachi and Robert L. Pego. Asymptotic stability of toda lattice solitons. Preprint, 
2007. 



